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ABSTRACT 

In  this  lecture,  we  present  some  suboptimum  STAP  algorithms  for  radar,  applied  to  moving  target 
indication.  In  addition,  we  show  some  MTI  results  obtained  with  the  multi-channel  airborne  experimental 
radar  AER-II  of  FGAN-FHR. 


1  INTRODUCTION 

In  the  previous  lecture  [1],  it  was  shown  that  if  the  interference  situation  (clutter,  jamming,  noise) 
is  known,  the  optimum  space-time  adaptive  filter  vector  w  is  given  by 

w  =  //Q_1a(u0,/0),  (1) 

where  a(u0,/0)  denotes  the  given  space-time  signal  vector, 


Q  =  JE,{(c  +  j  +  n)(c  +  j  +  n)//}  (2) 

is  the  space-time  clutter-plus-jamming-plus-noise  covariance  matrix,  and  p  is  a  normalisation  constant 
which  can  be  chosen  arbitrarily. 

The  large  size  of  the  space-time  covariance  matrix  precludes  the  use  of  optimum  space-time  adaptive 
processing  in  an  operational  radar  system,  for  it  requires  both  an  enormous  computational  effort  and  more 
training  data  than  is  usually  available.  Hence,  the  dimension  of  the  filtering  problem  must  be  reduced. 
The  resulting  “suboptimum  processors”  will  generally  outperform  optimum  STAP  if  the  covariance  matrix 
is  estimated  from  the  data,  since  they  require  less  computational  effort  and  fewer  training  data. 


2  ALGORITHMS 

2.1  The  Dimension  of  the  Clutter  Subspace 

In  the  lecture  on  principles  of  adaptive  array  processing  [5],  subspace  methods  were  introduced  as  a 
technique  for  reducing  the  amount  of  training  data  required  for  estimating  the  adaptive  weights. 
This  requires  knowledge  of  the  dimension  dim  CSS  of  the  clutter  subspace. 
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For  a  side  looking  linear  array,  Brennan  and  Staudaher  [6]  found  that 

dim  CSS  «  int{7V  +  (M  - l)/},1 

where  N  is  the  number  of  array  channels,  M  is  the  number  of  pulses,  and 

2 vpT 


(3) 


(4) 


is  the  ratio  between  the  normalised  clutter  Doppler  frequency  fDT  and  the  spatial  frequency 
with  vp  the  platform  velocity,  T  the  pulse  repetition  interval,  and  d  the  inter-element  distance. 

In  general,  equation  (3),  which  is  also  known  as  Brennan’s  mle,  only  gives  a  lower  bound  on  the  number 
of  clutter  eigenvalues.  It  does  not  hold  for  forward  looking  planar  arrays,  as  is  illustrated  by  Figure  1, 
which  shows  the  eigenvalues  of  an  estimated  space-time  clutter-plus-noise  covariance  matrix,  normalized 
to  the  noise  power. 


zi 


cos/?, 


Eigenvalues  of  the  estimated  space-time  covariance  matrix 


Figure  1 :  Eigenvalues  of  the  Space-Time  Covariance  Matrix  for  a  Forward  Looking  Planar  Array. 

2.2  The  Displaced  Phase  Centre  Antenna  (DPCA)  Technique 

Historically,  the  displaced  phase  centre  antenna  (DPCA)  technique  was  the  first  -  and  simplest  -  approach 
to  space-time  processing  of  clutter  echoes.  It  was  developed  for  implementation  in  RF  technology,  long 
time  before  the  advent  of  digital  signal  processing.  As  illustrated  in  Figure  2,  DPCA  is  based  on  a  side 
looking  antenna  arrangement  with  two  (or  more)  phase  centres  displaced  along  the  flight  axis.  The  aircraft 
motion  can  then  be  compensated  by  choosing  the  pulse  repetition  frequency  (PRF)  so  that  the  second 
phase  centre  assumes  the  position  of  the  first  phase  centre  after  one  pulse  repetition  interval  (PRI). 
Therefore,  any  two  successive  pulses  received  by  the  two  different  antenna  parts  appear  to  come  from  one 
phase  centre  fixed  in  space.  A  simple  two-pulse  canceller  subtracts  the  second  echo  received  by  the 
trailing  antenna  from  the  first  echo  received  by  the  leading  antenna. 


1  A  similar  expression  was  found  empirically  by  Klemm  [7]  almost  ten  years  earlier. 
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It  should  be  noted  that  the  DPCA  technique,  being  a  deterministic  space-time  processing  technique, 
requires  identical  channels,  antenna  phase  centres  displaced  exactly  along  the  flight  axis,  and  a  PRF 
perfectly  matched  to  the  platform  velocity  in  order  to  work  properly.  However,  due  e.g.  to  channel  errors 
and  crab  angle,  these  requirements  will  hardly  be  met  in  practice.  Adaptive  techniques  have  been  devised 
to  overcome  this  problem. 


Figure  2:  The  Displaced  Phase  Centre  Antenna  (DPCA)  Technique. 


2.3  Post-Doppler  STAP 

Post-Doppler  space-time  adaptive  processing  techniques  operate  in  the  range-Doppler  domain,  i.e.  after 
Doppler  filtering  has  first  been  applied.  The  primary  rationale  for  this  is  the  dependence  of  the  main  beam 
clutter  Doppler  frequency  on  the  cone  angle  (cf.  [1]),  which  causes  the  clutter  echoes  to  occupy  only  a  few 
Doppler  bins  (while  being  present  in  all  pulses),  as  is  illustrated  in  Figure  3.  This  Doppler  subbanding 
allows  to  use  smaller  adaptive  filters,  which  require  less  computational  effort  and  fewer  training  data. 
In  addition,  post-Doppler  methods  are  less  sensitive  to  element  errors. 


Clutter  power  in  the  sum  channel  [dB] 


Clutter  power  in  the  sum  channel  [dB] 


0.02  0.03  0.04 

Time  [s] 


Figure  3:  Clutter  Power  [dB]  in  the  Sum  Channel  of  a  Forward  Looking  Radar,  Plotted  Over 
Range  and  Pulse-to-Pulse  Time  Resp.  Doppler  Frequency  Modulo  PRF. 
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Post-Doppler  Adaptive  Beamforming 

The  simplest  post-Doppler  STAP  technique  is  post-Doppler  adaptive  beamforming ,  which  has  also  been 
called  frequency  dependent  spatial  processing,  or  factored  post-Doppler.  In  this  method,  the  NM  - 
dimensional  space-time  filtering  problem  is  divided  into  M  separate  N  -dimensional  adaptive 
beamforming  problems,  as  shown  in  Figure  4.  Both  the  computational  effort  and  the  amount  of  training 
data  required  are  thus  reduced  considerably. 


K  range  gates 

◄ - ► 


(separately  in  each  Doppler  bin) 

Figure  4:  Post-Doppler  Adaptive  Beamforming. 

Since  post-Doppler  adaptive  beamforming  is  only  spatially  adaptive,  the  Doppler  bins  must  be 
decorrelated  in  order  for  this  technique  to  work  well.  As  the  correlation  between  Doppler  bins  is  due  to 
mainbeam  clutter  leakage  into  the  Doppler  sidelobes,  it  can  be  reduced  by  using  a  windowed  Fourier 
transform,  which  leads  to  low  Doppler  sidelobes.  However,  the  heavy  tapering  required  also  causes  a 
broadening  of  the  Doppler  mainlobe  and  significant  taper  loss. 

Post-Doppler  adaptive  beamforming  is  suitable  for  moving  target  detection  using  a  multi-channel 
synthetic  aperture  radar,  for  the  long  pulse  trains  transmitted  lead  to  high  Doppler  resolution. 


Adjacent-Bin  Post-Doppler  STAP 

The  poor  performance  of  post-Doppler  adaptive  beamforming  can  be  improved  by  providing  the 
architecture  with  temporal  adaptivity.  One  way  this  can  be  achieved  is  by  adaptively  processing  the  spatial 
returns  from  several  adjacent  Doppler  bins  ( adjacent-bin  post-Doppler  STAP  [8]).  In  this  method,  which  is 
also  called  extended  factored  STAP,  the  NM  -dimensional  space-time  filtering  problem  is  divided  into 
M  separate  NL  -dimensional  space-time  filtering  problems,  as  shown  in  Figure  5.  Both  the 
computational  effort  and  the  amount  of  training  data  required  are  thus  reduced  considerably. 
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Adjacent-bin  post-Doppler  STAP,  while  performing  much  better  than  post-Doppler  adaptive 
beamforming,  is  sensitive  to  Doppler  tapering,  which  increases  the  dimension  of  the  clutter  subspace  [1]. 


PRI-Staggered  Post-Doppler  STAP 

PRI-staggered  post-Doppler  STAP  [6], [9]  uses  a  bank  of  Mf -pulse  Doppler  filters  that  produces 
L  —  M  —  M'+ 1  output  pulses  for  each  filter,  see  Figure  6.  The  spatial  returns  from  the  L  output  pulses 
are  then  processed  adaptively.  As  in  adjacent-bin  post-Doppler  STAP,  the  NM  -dimensional  space-time 
filtering  problem  is  again  divided  into  M  separate  NL  -dimensional  space-time  filtering  problems. 


M  pulses 


X" 


>-r 


K  range  gates 


Doppler  filter  bank 
(Mf  pulses) 

K  range  gates 


Space-time  adaptive 
processing 


Figure  6:  PRI-Staggered  Post-Doppler  STAP. 
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Given  the  same  number  of  degrees  of  freedom  and  Doppler  sidelobes,  the  performance  of  PRI-staggered 
post-Doppler  STAP  is  usually  better  than  that  of  adjacent-bin  post-Doppler  STAP  [1]. 


Symmetric  Auxiliary  Sensor  /  Echo  Processor 

A  space-time  adaptive  processor  closely  related  to  PRI-staggered  post-Doppler  STAP,  and  particularly 
well-suited  for  linear  arrays,  is  described  in  [10].  It  is  called  symmetric  auxiliary  sensor  /  echo  processor 
and  makes  use  of  the  fact  that  the  concepts  of  overlapping  uniform  subarrays  and  symmetric  auxiliary 
sensors,  which  are  illustrated  in  Figure  7,  lead  to  the  same  output  but  different  processing  schemes.  In  fact, 
the  only  difference  between  the  two  concepts  is  that  all  those  central  elements  common  to  both  subarrays 
are  pre-summed.  Analogously,  the  concepts  of  overlapping  time  intervals  (which  is  used  in  PRI-staggered 
post-Doppler  STAP)  and  symmetric  auxiliary  echoes  lead  to  the  same  results  but  different  processing 
schemes  (in  this  case,  the  sum  beamformers  in  Figure  7  must  be  viewed  as  Doppler  filters). 


Overlapping  time  intervals 


Symmetric  auxiliary  sensors 

(Summation  of  central  elements) 

Symmetric  auxiliary  echoes 


Figure  7:  From  Overlapping  Subarrays  to  Auxiliary  Sensors. 


As  shown  in  Figure  8,  the  symmetric  auxiliary  sensor  /echo  processor  [10]  uses  both  a  central 
beamformer  and  a  central  Doppler  filter,  along  with  auxiliary  sensors  and  auxiliary  echoes,  thus  reducing 
the  dimension  of  the  adaptive  problem  -  and  both  the  computational  effort  and  the  amount  of  training  data 
required. 


K  range  gates 

◄ - ► 


Space-time  adaptive 
processing 


Figure  8:  The  Symmetric  Auxiliary  Sensor  /  Echo  Processor  (ASEP). 

For  linear  arrays,  the  structure  of  the  symmetric  auxiliary  sensor  /  echo  processor  leads  to  an  efficient 
algorithm  for  determining  the  inverse  of  the  space-time  clutter-plus-noise  covariance  matrix. 
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Joint  Domain  Localized 

As  discussed  at  the  beginning  of  this  section,  the  rationale  behind  post-Doppler  STAP  is  the  fact  that  in 
view  of  the  dependence  of  the  main  beam  clutter  Doppler  frequency  on  the  cone  angle,  clutter  suppression 
in  the  range-Doppler  domain  is  less  complex  than  in  the  range-(pulse)  time  domain.  In  fact,  this  Doppler 
subbanding  allows  to  use  smaller  adaptive  filters,  which  require  less  computational  effort  and  fewer 
training  data. 


Similarly,  for  array  antennas  with  a  large  number  of  spatial  channels,  using  beams  formed  from  the 
outputs  of  these  channels  may  be  more  economical.  Such  techniques  are  called  beamspace  STAP. 


Joint  domain  localized  (JDL,  [11])  is  probably  the  best-known  beamspace  post-Doppler  STAP  technique. 
In  this  method,  the  NM  -dimensional  space-time  filtering  problem  is  divided  into  N'  M  separate 
L' L  -dimensional  space-time  filtering  problems,  as  shown  in  Figure  5.  Both  the  computational  effort  and 
the  amount  of  training  data  required  are  thus  reduced  considerably. 
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Figure  9:  Joint  Domain  Localized  (JDL). 


2.4  Data-Dependent  STAP 

The  rationale  behind  data-dependent  STAP  is  the  idea  that  adaptively  reducing  the  dimension  of  the 
filtering  problems  will  lead  to  a  more  efficient  use  of  the  available  degrees  of  freedom.  We  will  discuss 
these  techniques  only  very  briefly,  as  they  are  not  specific  of  STAP.  For  details,  the  reader  is  referred 
to  [3]. 


Principal  Components  Inverse  (PCI) 

Principal  components  inverse  (PCI)  and  its  modifications  are  subspace  methods  which  replace  the  inverse 
of  the  space-time  clutter-plus-noise  covariance  matrix  by  a  weighted  projection.  It  can  be  interpreted  as  a 
lean  matrix  inversion  (LMI)  method  [5]. 
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Multistage  Wiener  Filter 

The  multistage  Wiener  filter  (MWF,  [12])  selects  a  basis  of  the  clutter  subspace  by  successively  choosing 
rank-one  vectors  so  that  the  mean  square  error  is  minimized,  but  no  knowledge  of  the  full  space-time 
covariance  matrix  is  required. 

2.5  The  Space-Time  Adaptive  FIR  Filter 

Using  the  Cholesky  decomposition 


Q  1  =  R^R  (5) 

of  the  space-time  clutter-plus-noise  covariance  matrix,  the  output  of  the  optimum  adaptive  space-time 
filter  can  be  written  as 


w"z  =  (//•Q"1a(u0,/0)f  z 

=  //-a"(u0,/0)R"Rz  (6) 

=  //-(Ra(u0,/0)f  Rz. 

The  triangular  matrix  R  decorrelates  ( whitens )  the  received  signal  and  matches  the  space-time  steering 
vector  a(u0,/0)  .  This  whitening  resp.  matching  can  also  be  effected  using  a  finite  impulse  response  (FIR) 
filter. 

The  space-time  adaptive  FIR  filter  [13]  is  an  extension  of  the  linear  prediction  error  filter  to  space-time 
signals.  This  filter  can  be  used  for  minimizing  the  interference  power  in  the  received  signal  because  it 
minimizes  the  prediction  error  in  the  least  squares  sense. 

The  coefficients  of  a  space-time  FIR  filter  operating  in  the  temporal  dimension  are  given  by  the  first  block 
column  of  the  inverse  of  the  covariance  matrix  QL  of  the  space-time  data  segment  of  temporal  length  L  . 
By  including  the  spatial  matching  (secondary  spatial  beamforming)  in  the  filtering  operation,  one  obtains  a 
filter  vector  of  length  NL ,  where  N  is  the  number  of  spatial  channels.2 

The  space-time  FIR  filter  offers  several  advantages: 

•  As  the  temporal  length  L  of  the  filter  can  be  chosen  independently  of  the  system  dimensions 
(number  of  channels  N  ,  number  of  pulses  M ),  the  number  of  filter  coefficients  can  be  small,  so 
that  few  training  data  are  required,  and  the  filtering  operation  can  be  carried  out  in  real-time. 

•  The  filter  can  also  be  used  with  inhomogeneous  clutter  or  with  a  manoeuvring  platform;  however, 
it  will  then  have  to  be  re-adapted  from  pulse  to  pulse. 

On  the  downside,  the  filtering  operation  causes  some  signal  loss,  as  it  shortens  the  data  sequence  from 
M  pulses  to  M  —  L  + 1  pulses.  This  effect  is  visible  in  Figure  10,  which  shows  the  performance  of  the 
space-time  adaptive  FIR  filter  (< dashed  line )  for  a  forward  looking  planar  array  antenna.  Compared  to 
optimum  STAP  using  the  subarray  outputs  {solid  line),  the  space-time  adaptive  FIR  filter  shows 
some  losses  in  pass  band,  but  achieves  the  same  narrow  clutter  notch.  Conventional  processing, 
i.e.  beamforming  plus  Doppler  filtering  (< dotted  line),  shows  significant  losses  in  pass  band  and  does  not 
achieve  a  narrow  clutter  notch. 


2  For  details,  the  reader  is  referred  to  [2]. 
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Performance  of  the  space-time  adaptive  FIR  filter 


Figure  10:  The  Performance  of  the  Space-Time  Adaptive  FIR  filter  (solid  line) 
Compared  to  Optimum  STAP  Using  the  Subarray  Outputs  (dashed  line) 
and  Beamforming  +  Doppler  Filtering  (dotted  line). 


3  EXPERIMENTAL  RESULTS 

The  results  shown  in  this  chapter  were  obtained  using  data  acquired  with  the  Airborne  Experimental 
Radar  (AER-II)  of  FGAN-FHR.  This  quasi-operational  multi-channel  X-band  SAR/MTI  radar  system 
features  an  active  electronically  scanned  array  antenna.  It  has  now  been  succeeded  by  the  Phased  Array 
Multifunctional  Imaging  Radar  (PAMIR). 

Figure  11  shows  the  outcome  of  a  controlled  SAR/MTI  experiment  with  eight  slowly  moving  ground 
targets.  In  both  subplots,  the  signal  power  (in  dB)  is  plotted  over  the  ground  coordinates.  While  the  targets 
are  hidden  by  clutter  in  the  sum  beam  (left  subplot),  they  are  clearly  visible  after  clutter  suppression  using 
post-Doppler  adaptive  beamforming  (right  subplot). 


Figure  11:  The  Result  of  Space-Time  Adaptive  Processing: 
Signal  Power  [dB]  Plotted  Over  (x,y)  (Image  courtesy  J.  Ender). 
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Figure  12  shows  the  result  of  combining  a  SAR  image  of  the  scene,  in  which  the  moving  targets  are  not 
visible,  with  the  GMTI  information  shown  in  the  previous  image. 


Figure  12:  Fusion  of  SAR  Image  and  GMTI  Information  (Image  courtesy  J.  Ender). 


Figure  13  shows  the  results  of  another  SAR/MTI  experiment,  using  targets  of  opportunity.  Cars  on  a 
motorway  junction  were  detected  and  correctly  positioned;  for  some  of  them,  tracks  could  be  established. 


Figure  13:  SAR  Image  of  a  Motorway  Junction,  Overlaid  with  the  Accumulated  Moving  Target 
Detections  (left)  and  those  Targets  for  which  Tracks  could  be  Established  (right) 

(Image  courtesy  J.  Ender). 
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